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We derive a lower bound for the concurrence of mixed bipartite quantum states, valid in arbitrary 
dimensions. As a corollary, a weaker, purely algebraic estimate is found, which detects mixed 
entangled states with positive partial transpose. 
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In classical physics, one can always divide a system into 
subsystems, such that complete information on the entity 
implies a complete description of its individual parts, and 
vice versa. In quantum physics, this no longer holds true: 
whilst one can still divide a system into subsystems, a 
complete description of the system state in terms of a 
pure state does not necessarily assign a pure state to 
each subsystem. The subsystems of generic pure states 
are correlated in a way without classical analog - they 
are entangled. 

While such quantum entanglement arguably incar- 
nates the key difference between the quantum and the 
classical world, and is nowadays understood as a resource 
in various tasks of quantum information processing 
such as cryptography, teleportation, and quantum com- 
putation, it remains hard to quantify, for arbitrary quan- 
tum states In particular, when coupled to an envi- 
ronment, pure quantum states rapidly evolve into mixed 
states which bear entanglement together with classical 
probabilistic correlations, and the latter have to be dis- 
tinguished from the former. Furthermore, the complete 
characterization of the nonclassical correlations of a given 
state becomes an ever more complex task as the Hilbert 
space dimension increases, thus turning into a computa- 
tionally extremely intricate problem. 

No equally versatile as computationally manageable 
entanglement measure for mixed states is available so far, 
although various more or less pragmatically motivated 
quantities have been proposed. The most popular indica- 
tor of entanglement is the positive partial transpose (ppt) 
criterion Q and variants thereof, such as negativity 
though these do not reliably detect arbitrary entangled 
states. Another approach for quantifying entanglement is 
through entanglement witnesses Q which, however, need 
to be constructed anew for each given quantum state, 
and such construction can be rather involved [^. Fi- 
nally, there are mixed state generalizations of pure state 

entanglement measures Hi ma ins, which, m gen- 
eral, require a high dimensional optimisation procedure. 
By construction, any numerical evaluation of these latter 
quantities only yields upper bounds for the entanglement 
of a given state but cannot reliably distinguish it from 
separable states, let alone provide a reliable quantitative 
estimate of the state's actual degree of entanglement. 

In the present Letter, we improve on that situation: 
We derive a /ower bound of concurrence ~ ^ quantity 



which is strictly larger than zero for nonvanishing entan- 
glement - of mixed bipartite quantum states in arbitrary 
dimensions. Our bound is given by a purely algebraic ex- 
pression which is readily evaluated for arbitrary states, 
and can be tightened numerically on a relatively low- 
dimensional parameter space, of reduced dimension as 
compared to hitherto available optimisation procedures 
h3. This complements already available upper bounds 
|lOl [ll| and provides, for the first time, a rather precise 
estimate of the actual value of concurrence. 

We start out with the definition of a pure state's 
concurrence as c{'4)) = ~ Trg^, where the re- 

duced density matrix is obtained by tracing over one 
subsystem. The concurrence of mixed states g is then 
given as the convex roof 

c{g)^mi^p,c{-^.,), g^^p,\^,){^,\, > 0, (1) 



of all possible decompositions into pure states Con- 
sequently, c{g) vanishes if and only if g exhibits purely 
classical correlations, i.e. if the state is separable and 
hence can be represented as a convex sum over product 



states, g — J^iPi^'f'' ® £'^^^ '^i*^ Pi ^ Oj ^^^^ g\ 

(2) 

g] states on the subcomponents Tii and 7^2 of the total 
Hilbert space TL = Tii ® 0,2- Given the dimensions ni 
and n2 of Hi and H2, respectively, eq. (1) defines a high 
dimensional optimisation problem which is rather cum- 
bersome to solve. Furthermore, as already mentioned 
above, such optimisation can only yield an upper bound 
for c(p) , by virtue of the definition of the infimum. 

To estimate c{p) from below, we first replace, for con- 
venience, the l^'j) by the subnormalized states \ipi) — 
y/pi\'^i) in eq. (0). Given a vahd decomposition of 
g into subnormalized states, any other suitable set 
is obtained ^3 by transformations V € (jNxr^ with r 
and N the lengths of the sets and {j-f/'i)}, respec- 

tively. 



J2) 



and 



N 



(2) 



It is now crucial to realize that the concurrence of a pure 
state \ip) can be expressed as the square root of the func- 
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tion 



/('01,l/'2,l/'3, V'4) = (V'2|V'l)(V'4|V'3) - 

Tri((Tr2|V'i)(^2|)(Tr2|V'3)(V'4|)) 



(3) 



evaluated at = V'l = V'2 = V'a = V'47 where Tri and Tr2 
denote the traces over the first and the second subsystem. 
/ is linear in its first and third, and anti-linear in the 
second and fourth argument. Due to these properties 
the definition |^ can be reformulated as an infimum over 
transformations V: 

N 1 

c{g) = MC, with C ^^(^[V(E)VAV^ (E)V'']'^y . 

1=1 

(4) 

Herein, the tensor A, defined by A'™ = f{<j)j,(j)i,(j>k,(j>m) 
[T^ is hermitian, A'^JJ} = {Aj^)*, and symmetric with 
respect to a simultaneous exchange of both its co- and 
contravariant indices Ap! = A^j. Due to the symme- 
try of the transformation V (E) V under exchange of the 
subsystems of A, we can replace A^JJ! in eq. Q by the 
symmetrised elements 



1 



- ( 4'™ 



47) 



(5) 



which is equivalent to a symmetrisation over both sub- 
systems in eq. It can be shown that A is pos- 
itive semidefinite and that its support lies in an anti- 
symmetric subspace, i.e., all elements of A with respect 
to fully symmetric linear combinations of product states 
vanish. Since the antisymmetric subspace has dimension 
m — ni{ni — l)rt2('T-2 ^ l)/4, A has at most m non- 
vanishing eigenvalues. 

Due to the discussed symmetries A can be expanded 
in a basis of real symmetric matrices A" G pj^'X'" 



4" = E^"/3Aifc< 



(6) 



with B hermitian and positive semi definite. With the 
eigenvalues and associated eigenvectors of B {Bx°' — 
^aX°' and — [xf, ... ,xf ,...]) we can construct a prop- 
erly normalised eigensystcm T" of A 

To. = ^Y.^P'''^"^^ = ^"6"^° , a = 1, . . . , m. (7) 



We explicitly take into account the free phase factors 
exp(i(^Q,), as they will be crucial in the following. Conse- 
quently, 



j\lni \ ^ rpa frpa \* \ ^ /-ra ('~rCL 

•^jk ^ 2^ ^jk K-'-lrn) ~ ^ jk [■'Im 

a a 

Hence, eq. Q can now be rewritten as 



(8) 



(9) 



the infimum of which gives the concurrence of the mixed 
state g. 

Note that © resembles the concurrence vector in- 
troduced in a quantity with elements analogous to 

Ez=i In O it was shown that the con- 

currence vector vanishes identically, for suitably chosen 
V, if and only if g is separable, thus providing a sep- 
arability criterion. Whilst the very same equivalence 
holds for vanishing C, the explicit expression allows 

us to proceed further: The Cauchy-Schwarz inequal- 
i_ 1 

ity (EaS^a)' > Ea^^aya, and Eal^^l > 

I Y,a ^"1' s imply 



N 



i V 

c{g)>Mj2\[v(j2zo,T-] 

i—1 a 



(10) 



for any set Za = yaexp{i(l)a), with ?/„ > 0, J^aUl = 1- 
The infimum on the rhs is given by ^i—^i^i K, where Xj 
are the singular values of T = z^T^, i.e., the square 
roots of the eigenvalues of the positive hermitian matrix 
TT^ in decreasing order |lOj |. Hence, we arrive at the 
desired lower bound, 



c{g) > Ai -^Ai, 



(11) 



with the Xj dependent on the choice of the and 0q. 

Note that each set {ya,4'a} provides a lower bound 
of c(p), which can be tightened by numerical optimiza- 
tion. However, all the examples we have considered so 
far suggest that there is one matrix T" that gives the 
main contribution to the rhs of eq. H1U|I . Hence, the 
singular values of this matrix provide a purely algebraic 
lower bound for c, which often leads to satisfactory re- 
sults even without further numerical refinement. As an 
example, consider the family of two spin-1 states 



1 



Qa = 



l + 8a 



a e [0, 1] 



with 13 = (l + a)/2 and 7 = y/l ~ 



(12) 

described in [l6l| : 
The state ga has positive partial transpose in the entire 
range of a, i.e., the standard criterium ^ identifying a 
mixed state as nonseparable via negativity of its partial 
transpose is inoperational here. Notwithstanding, the 
singular values of associated with the largest eigen- 
value of B in eq. © above already provide a positive 



lower bound 



.(1) 



= -^i^^ - J2i>i K''' for c{ga), as shown 
in Fig. n ga is detected as entangled without any need 
for further numerical optimisation. However, the figure 
also shows that taking into account all T", numerical 
optimisation significantly raises the lower bound. The 
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FIG. 1: Lower bound (full line) of the concurrence c of the 
family of two spin-1 states Qa Ib], together with its purely al- 
gebraic approximation cp' = A^^' — Xl^^j^ A^^' (dashed line), 
over the entire parameter range a = ... 1. Whilst Qa has 
positive partial transpose and is therefore not detected via 
the negative partial transpose criterion for entanglement, al- 
ready the algebraic approximation of our lower bound inden- 
tifies the state as nonseparable, without need for numerical 
optimisation, for all a. Optimal lower bound and optimal up- 
per bound (dotted line, obtained by minimizing C in ||1J 
confine the actual value of concurrence to an interval with 
relative uncertainty of the order of approx. 10%. 

remaining gap with respect to the upper bound obtained 
from numerical minimization of C in Q leaves a relative 
uncertainty of the order of approx. 10% on the actual 
value of c{q), in the particularly pathological test case 
considered here. 

Let us finally note that the original definition of con- 
currence (9,) is nicely embedded in our formalism. In the 
case of two-level systems one has m = 1, i.e., there is 
only one non- vanishing matrix T^. Therefore, eq. 



simplifies to 

N 

C = Y^\[VT^V-U , (13) 

i=l 

and the infimum can be derived analytically. Indeed, it 
is found that coincides with r defined in the original 
work 10]. 

In conclusion, we have shown that a suitable repre- 
sentation of the concurrence of bipartite mixed quan- 
tum states in terms of the eigensystem of a tensorial 
quantity allows for the derivation of a lower bound of 
c(p), for arbitrary p. Not only can this bound be tight- 
ened by an optimisation under the comparatively simple 
constraint = Ij over a complex vector space 

of dimension niirii — I)n2(n2 — l)/4 - at least a fac- 
tor 4nin2 smaller than dimensions of optimisation pro- 
cedures hitherto available. It also can be reduced to a 
purely algebraic bound which appears to provide good 
estimates, according to numerical data which comple- 
ment our analytical work. Last but not least, our re- 
sult can serve to derive lower bounds on the entangle- 
ment of formation i? Q of a mixed quantum state, thus 
quantifying the minimally required resources to prepare 
q: Given any monotonously increasing, convex function 
£ which satisfies £{c{%}])) < —Trgrloggr, it follows that 
E{g) > £{c{g)), with the rhs bounded from below by our 
bound ((TT|l . 
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